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ABSTRACT 


Thermodynamic  and  transport  properties  of  rubidium  andcesium 
vapors  are  determined  from  theory,  supplemented  by  experimen¬ 
tal  results.  The  sitxiation  for  which  results  are  obtained  involves 
negligible  ionization  of  the  alkali  atoms. 

The  model  used  is  an  imperfect  gas  obeying  the  virial  equation  of 
state.  The  second  virial  coefficient  is  considered  adequate  to  ac¬ 
count  for  deviations  from  the  ideal  behavior  in  the  temperature 
and  pressure  range  under  consideration.  It  is  assumed  that  the 
vapors  are  comprised  of  alkali  atoms  whose  interactions  are  tne 
pairwise  Lennard- Jones  (6-12)  type. 

On  the  basis  of  the  two  parameters  which  occur  in  the  Lennard- 
Jones  potential  for  interactions  between  some  alkalis  and  noble 
gae  atoms,  and  the  known  parameters  for  the  noble  gas  atoms, 
the  parameters  are  estimated  for  the  alkalis.  From  these  pa¬ 
rameters  the  desired  thermodynamic  and  transport  properties  of 
rubidium  and  cesium  vapors  are  obtained.  Numerical  values  are 
given  for  the  important  properties. 


I. 


INTERMOLECULAR  POTENTIALS  FROM  SPECTRAL  SHIFTS 


Foreign  gases  can  cause  a  shift,  commonly  called  the  pressure  shift,  in 
the  central  frequency  of  the  spectral  distribution  from  radiating  atoms.  In  the 
excited  state,  the  radiating  atom  possesses  a  different  energy  of  interaction 
with  another  atom  than  when  in  the  ground  state.  This  difference  in  energy 
can  give  rise  to  a  frequency  shift  and  broadening  of  the  spectral  line.  In  the 
case  of  some  foreign  gases  the  average  mutual  potential  energy  can  be  greater 
when  the  radiating  atom  is  in  the  excited  state;  for  other  foreign  gases  the 
reverse  is  true.  Hence,  the  frequency  shift  may  be  positive  or  negative. 
Measured  values  of  this  frequency  shift  provide  information  regarding  inter- 
molecular  potentials. 

A  simple  model  of  the  origin  of  the  pressure  shift  is  shown  in  Eqs.  (1) 
through  (6).  The  frequency  emitted  by  a  radiating  atom  unperturbed  by 
foreign  gases  is 


V 

o 


(1) 


where  (E^  -  E^)  is  the  energy  difference  between  the  excited  and  ground  states 
which  gives  rise  to  the  frequency  v^.  In  the  presence  of  a  foreign  atom,  both 
excited  and  lower  levels  of  the  emitting  atom  will  be  perturbed,  so  that  the  new 
emitted  frequency  is 
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h  ■  h 


(2) 


where  V2(r)  is  the  interatomic  interaction  for  the  excited  state  and  V^(r) 
refers  to  the  ground  state.  The  preceding  equation  may  be  written  as 


V  ^  =  Av 

o 


V^(r)  -  V^(r) 
h 


(3) 
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where  Av  is  the  frequency  shift  caused  by  the  foreign  atoms.  If  there  are  m 
foreign  atoms  at  a  distance  r^  from  the  radiating  atom,  the  Lennard- Jones 
potential  gives 
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(5) 


In  Eqs.  (4)  and  (5),  p(r)  is  the  probability  that  the  radiating  atom  and  the 
foreign  atom  are  a  distance  r  apart.  In  general,  p(r)  will  have  the  following 
form: 


p(r) 


A.-V(r)/kT 


(6) 


where  A  is  a  normalization  constant. 

From  the  measured  values  of  frequency  shift  in  the  hyperfine  spectra  of 
alkalis,  as  well  as  the  temperature  coefficients  of  the  shift,  this  writer^  has 
obtained  the  parameters  in  the  Lennard- Jones  potential  for  the  following 
interactions:  rubidium-helium,  rubidium- neon,  cesium-helium,  and  cesium- 
neon.  The  temperature  coefficient  of  the  shift  is  (Av)-  Using  a  more 
detailed  theory  than  that  suggested  by  Eq.  (4),  the  results  given  in  Table  I 
were  obtained. 


2 


IL  LENNARD-JONES  (6-12)  PARAMETERS  FOR  RUBIDIUM 
VAPOR  AND  CESIUM  VAPOR 

Very  little  information  exists  regarding  the  form  of  an  intermolecular 
potential  which  exists  between  alkali  atoms.  Because  of  the  increasing  impor¬ 
tance  of  alkalis  in  modern  technology,  it  is  desirable  to  have  some  knowledge 
of  this  intermolecular  potential.  The  Lennard- Jones  (6-12)  potential  is 
reasonably  satisfactory  in  many  cases.  However,  in  the  interpretation  of  some 
phenomena,  including  alkalis  and  large  rare  gas  atoms,  the  writer  has  found 
this  potential  unsatisfactory  in  a  quantitative  way-  No  conclusive  evidence  has 
been  found  by  the  writer  of  the  usefulness  of  the  Lennard- Jones  potential  in 
alkali-alkali  interactions.  In  this  paper,  it  will  be  accepted  as  a  satisfactory 
app  roximation . 

Some  published  data  regarding  the  effects  of  foreign  gases  on  the  spectral 
lines  of  rubidium  and  cesium  can  be  reinterpreted  to  obtain  the  constants  in 
the  Liennard- Jones  potential  for  these  alkalis.  The  published  results  can  now 
be  looked  upon  in  a  slightly  different  fashion  to  obtain  the  Eennard- Jones  pa¬ 
rameters  for  rubidium- rubidium  and  cesium-cesium,  as  explained  below. 

The  L»ennard- Jones  potential  for  the  interaction  of  two  different  types 
of  molecular  species,  i  and  j,  separated  a  distance  r,  is  given  by 


V.. 

ij 


(7) 


where  all  symbols  have  their  usual  significance.  The  attractive  (negative) 
term  is  the  London  dispersion  term  in  the  van  der  Waals  energy.  This  dis¬ 
persion  energy  term  can  be  written  more  compactly  as 


V.. 
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(attraction) 


(8) 
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The  various  approximate  forms  which  have  been  used  for  the  C. .  give  rise  to 
2 

the  inequality  that 


<  C.. 


—  11 


C.. 

JJ 


(9) 


or  equivalently 


c 


ij 


(e..  €..) 
11  ij 


1/2 


(10) 


The  subscripts  ii  and  jj  refer  to  interactions  of  like  pairs  of  molecules.  One 
would  expect  the  equality  in  Eqs.  (9)  and  (10)  to  be  approached  more  closely 
as  the  molecules  become  more  and  more  alike.  Alkali  and  noble  gas  atoms 
are  as  unlike  as  atoms  can  be,  hence  the  equality  in  Eqs.  (9)  and  (10)  is  not 
to  be  expected. 

equation  (10)  is  commonly  used  as  a  basis  for  assigning  parameters  for 
unlike  interacting  partners  when  the  parameters  for  the  individual  species  are 
known  separately.  There  seems  to  be  no  theoretical  basis  for  the  usual  assign¬ 
ments  that 


and  that 


€..  =  (e..  €..) 

ij  11  jj' 


1/2 


<r. .  =  rr  (<r  .  +  <r. .) 
ij  2  '  11  jj' 


(11) 


(12) 
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Equations  (11)  and  (12)  are  reasonably  satisfactory  for  sinailar  substances. 

In  this  case  Eq.  (11)  will  be  used*  but  at  the  same  time  the  inequality  in  Eq.  (10) 

will  be  maintained.  Hence  only  a  lower  bound  on  <r..  can  be  assigned  if  both 

J  J 

<r. .  and  <r..  are  known.  The  results  are  far  more  sensitive  to  the  <r's  than 
ij  11 

to  the  c's.  Therefore  it  seems  probable  that  discrepancies  between  theory 
and  experiment  are  the  result  of  uncertainties  in  the  (r's  rather  than  in  the  c 's. 
When  the  equalities  in  Eqs.  (11)  and  (12)  fail,  variations  in  c's  are  usually 
not  large  enough  to  cause  sizeable  discrepancies.  Hence  the  choice  of  Eqs  (10) 
and  (11),  rather  than  (11)  and  (12),  as  the  basis  of  analysis  seems  reasonable. 

Table  I  shows  the  results  of  the  above  calculations.  Eqiiation  (11)  is 
used  to  determine  the  c's  for  rubidium- rubidium  and  cesium-cesium  interactions. 
If  both  Eqs.  (10)  and  (11)  are  used,  the  following  relationship  is  found  as  a 
lower  bound  for  o*..: 


<r. . 


JJ 


(13) 


The  parameters  obtained  in  Table  I  can  be  checked  by  calculating  the  London 
dispersion  parameters  and  comparing  them  with  the  recent  calculations  of 

3 

P.  R.  Fontana  for  both  Rb-Rb  and  Cs-Cs.  The  results  of  Table  I  yield  for 
Rb-Rb 


4ct  =  673  X  10 


.60 


ere  cm 


and  for  Cs-Cs 


4c  <r 


6 


807  X  I0"^° 


6 

erg  cm 
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Table  I.  Parameters  in  the  Lennard- Jones  (6-12)  Potential 


Interaction 

O’  X  10®  cm 

(  X  10^  ev 

Reference 

He  -  He 

2.  56 

0. 0882 

Footnote  2 

Ne  -  Ne 

2.  78 

0. 1406 

Footnote  3 

Rb  -  He 

3.  24 

0.416 

Footnote  1 

Rb  -  Ne 

3.43 

0.  536 

Footnote  1 

Rb  -  Rb 

4.  10 

1. 962 

From  Rb  -  He 

Rb  -  Rb 

4.  23 

2.  031 

From  Rb  -  Ne 

Rb  -  Rb 

4.  17  ±  0.  07 

2.  00  ±  0.  03 

Average 

Cs  -  He 

3.  39 

0.  352 

Footnote  1 

Cs  -  Ne 

3.  58 

0.  458 

Footnote  1 

Os  —  Os 

4.49 

1. 406 

From  Cs  -  He 

Os  -  Os 

4.  61 

1. 450 

From  Cs  -  Ne 

Os  —  Os 

4.  55  ±  0.  06 

1. 423J:0.  02 

Average 

6 


The  above  calculations  are  based  on  the  average  values  of  c  and  o*.  Fontana's 
calculations  yield  for  Rb-Rb  and  Cs-Cs,  respectively 


A  6 

4c  S’  = 


1570  X  10 


-60 


ergs  cm 


and 


4c o’^  =  ZZOO  X  lO"^®  ergs  cm^ 

The  latter  values  cannot  give  both  <  and  o-  separately,  but  only  the  product 
as  indicated.  T£  one  assumes  that  the  values  of  c  are  correct  as  given  in  Table  I, 
Fontana's  result  would  yield 


<r  (Rb-Kb)  =  4.  8  x  10‘®  cm 


and 


<r  (Cs-Cs)  =  5.  3  X  10"®  cm 

In  both  cases  the  results  are  about  15  per  cent  larger  than  the  corresponding 
values  in  Table  I. 
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III.  THERMODYNAMIC  PROPERTIES 


Thermodynamic  properties  of  rubidium  and  cesium  vapors  are  readily 
obtained  by  considering  vapors  as  imperfect  gases  obeying  the  virial  equation 


of  state- 
energy), 


The  properties  below  are  given  per  mole:  S  (entropy),  E 


(heat  capacity  at  constant  pressure),  and  H  (enthalpy). 


(internal 


'N/ 


(S  -  S'')/P  +  {R/P)ln  (P/P”)  =  - 


dB 

dT 


(14) 


(E  -  E°)/P  =  -  T 


dB 

dT 


(15) 


(£p  -  = 


2 

d^B 


jdT^, 


(16) 


(H  -  H”)/P  = 


B  -  T 


dB 

dT 


(17) 


The  equations  above  apply  to  the  regions  of  temperature  and  pressure  where 
the  second  virial  coefficient  B  is  sufficient  to  represent  the  deviations  of  the 
alkali  vapors  from  ideal  behavior. 

Hirshfelder,  Curtiss,  and  Bird  have  calculated  B  ,  T  dB  /dT  ,  and 
^22^  ^2  ^ 

(T  )  d  B  /d(T  )  where  B  is  the  second  virial  coefficient  in  terms  of  the 
universal  parameter  T  ,  assuming  the  validity  of  the  Lennard- Jones  (6-12) 
potential.  The  parameter  T  =  kT/c  ,  where  €  is  defined  in  Eqs.  (7)  and  (11). 
The  actual  second  virial  coefficient  B  is  related  to  B  as  follows: 

B(T)  =  2/3  IT  N(r^B*^(T*),  where  <r  is  defined  in  Eqs.  (7)  and  (13)  and  N  is 
Avogadro's  number.  In  this  report,  the  quantities  appearing  on  the  right  hand 
side  of  Eqs.  (14),  (15),  (16),  and  (17)  are  listed  in  Table  II  for  rubidium  and 
cesium.  These  results  are  not  appreciably  different  from  those  reported  by 
Evans,  Jacobson,  Munson  and  Wagman,  which  were  obtained  in  a  different 
manner 
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Table  11.  Thermodynamic  Properties  of  Alkali  Vapors 


IV.  TRANSPORT  PROPERTIES 


Expressions  for  the  transport  properties  of  vapors  can  be  given  in  terms 

of  the  well-known  collision  integrals  (or  omega  integrals)  of  Chapman  and 

4  5 

Cowling  and  of  Hirshfelder,  Curtiss,  and  Bird.  The  expressions  for  the 

diffusion  coefficient  D,  the  viscosity  coefficient  ti,  and  the  heat  conduction 

coefficient  te  are  given  below  for  pure  substances.  These  expressions  are 

based  on  the  assumption  of  a  Lennard- Jones  (6-12)  interaction  between  molecules. 


D  = 


3 

IT 


1/2 


(2vk^T^/n) 
Ptto-^ 


(18) 


In  Eq.  (18),  |i  is  the  reduced  mass,  P  is  the  pressure,  and  <r  is  defined  in 

(s  t)* 

Eqs.  (7)  and  (13).  All  of  the  omega  integrals,  Q'  ’  ,  will  be  given  later. 


(irmkT)^^^ 

..^2,2)* 


In  Eq.  (19)f  m  is  the  mass. 


^  25  (irmkT)^^^  ^ 
32  ^  T] 


v.y2.2) 


(19) 


(20) 


In  Eq.  (20,  s  is  the  specific  heat  of  the  vapor  at  constant  volume.  The 
^  (s  t)* 

collision  integrals  ’  are  given  by 


12 


(8,t) 


* 


(21) 
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where  "i  =  (|i./2kT)^^^  is  the  reduced  (dimensionless)  initial  relative  speed  of 

the  colliding  partners,  with  actual  initial  relative  speed  v^. 

=  Zv  /  I(v  ,  e)  (1  -  cos®  6)  sin  0  dS  ,  (22) 

Jo  ^ 

In  Eq.  (22),  I(v^,6)  is  the  differential  scattering  cross  section.  The  collision 
integrals  n(s,t)*  can  be  evaluated  in  terms  of  a  single  \iniversal  parameter 
T  =  kT/c,  where  c  is  defined  as  in  Eqs.  (7)  and  (11).  Hirshfelder,  Curtiss, 
and  Bird  have  also  tabulated  the  pertinent  On  the  basis  of  these 

tabulations,  the  transport  coefficients  given  in  Eqs.  (18),  (19)i  and  (20)  can 
be  calculated  for  both  rubidium  and  cesium. 

In  this  paper  the  viscosity  coefficients  of  both  rubidium  vapor  and  cesium 

vapor  have  been  calculated  and  compared  with  one  other  known  calculation  for 

cesium.  In  a  paper  called  Some  Selected  Properties  of  Cesium,  ^  Irving  Granet 

calculated  transport  properties  of  cesium  vapor  in  a  straightforward  fashion. 

He  used  the  mean  free  path  theory  from  the  simple  kinetic  theory  of  gases. 

The  diameter  of  cesium  atoms  was  estimated  from  the  density  of  solid  cesium. 

Granet  assumed  that  cesium  vapor  is  essentially  monatomic.  This  is  certainly 

suggested  by  every  estimate  obtained  for  the  value  of  <r.  Granet  found  for 

cesium  ti  =  0.98x10  ^  sPt  poise.  The  present  calculations  give 

=  1 . 45  X  10  ^  {sIT/0^^*  ^^**')  poise  for  <r  =  4.  55  x  10  ®  cm  and 

t]  -  1.  12  X  10  ^  (•Tt/O^^*  ^^*)  poise  for  <r  =  5.  3  x  10  ®  cm.  The  collision 
(2  2)^ 

integral  R'  ’  '  is  of  the  order  of  unity.  The  present  calculation  agrees  quite 
well  with  the  results  of  Granet' s  simple  calculation. 

Table  III  lists  the  transport  properties  given  in  Eqs.  (18),  (19),  and  (20) 

4c  4c 

for  various  values  of  the  reduced  temperature,  T  .  To  convert  T  to  degrees 
Kelvin,  one  must  multiply  by  232  for  rubidium  and  by  165.  0  for  cesium.  A 
simplifying  assumption,  s^  -  3R/(2M),  has  been  made  in  Eq.  (20);  R  is  the 
molar  gas  constant  and  M  is  the  atomic  mass  number  (atomic  weight)  in  grams. 
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C  e  sium  Rubidium 


Table  III.  Transport  Properties  of  Alkali  Vapors 


-T,* 

T 

rt  X  10^ 

(poise) 

a  X  10® 

a 

PD 

b 

crl 

<r2 

al 

<tZ 

(tI 

<r2 

1 

1.  362 

1.028 

1.  187 

0.896 

4.01  X  10’^ 

3.03  X  10"^ 

1.  5 

2.01 

1. 520 

1.  756 

1.  325 

8.85  X  10"^ 

4.  39  X  10'^ 

2 

2.  60 

1.963 

2.  27 

1. 711 

1. 518  X  10"^ 

8.42  X  10"^ 

1  2.  5 

3.  13 

2.  36 

2.  73 

2.06 

2.28  X  10’^ 

1.  360  X  10’^ 

1  3 

3.  60 

2.  72 

3.  14 

2.  37 

3.  28  X  10’^ 

1. 981  X  10“^ 

1  3.5 

4.  04 

3.05 

3.  52 

2.66 

4. 14  X  10-1 

2. 70  X  10"^ 

4 

4.46 

3.  36 

3.  88 

2.  93 

5.  22  X  10"^ 

3.  51  X  10’^ 

5 

5.  21 

3.93 

4.  54 

3.43 

7.66  X  10“^ 

5. 38  X  10'^ 

10 

8.29 

6.26 

7.  23 

5.45 

2.46 

1.  943 

20 

13.  00 

9.81 

11.  34 

8.  56 

7.  77 

6.  81 

50 

23.  5 

17.  73 

20.48 

15.46 

35.43 

35.  50 

II 

b 

X  10  cm. 

=  4.8 

in-8 

X  10  cm 

1 

1.203 

0.  886 

0.675 

0.497 

1.621  X 10"^ 

1. 195  X  10"^ 

1.  5 

1.  779 

1.  311 

0.998 

0.  735 

3.  58  X  10"^ 

2.64  X  10"^ 

2 

2.  30 

1. 693 

1.288 

0.  950 

6.  14  X  10"^ 

4.  52  X  10"^ 

2.  5 

2.  76 

2.  03 

1.  549 

1.  141 

9.22  X  10"^ 

6.  80  X  10'^ 

3 

3.  18 

2.  35 

1.  785 

1.  315 

1. 277  X  10"^ 

9.41  X  10'^ 

1  3.  5 

3.  57 

2.63 

2.  00 

1.476 

1. 675  X  10’^ 

l.ZH  X  10'^ 

•  r4 

S  4 

3.  94 

2.  90 

2.  21 

1. 627 

2. 11  X  10'^ 

1. 557  X  10'^ 

^  5 

4.60 

3.  39 

2.  58 

1. 903 

3.  10  X  10'^ 

2.  28  X  10'^ 

10 

7.  32 

5.40 

4.  11 

3.03 

9.  94  X  10'^ 

7.  32  X  10"^ 

20 

11.49 

8.46 

6.44 

4.  75 

2.  93 

2.  16 

50 

20.  8 

15.  29 

11.64 

8  58 

14.  33 

10.  56 

ffj  =  4.  55 

in-8 

X  10  cm, 

=  5.  3 

10-8 

X  10  cm 

a  Units  are  cal  g  ^  sec  ^  deg  ^  K 
b  Units  are  cm^  sec”^  atm 
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Figure  1  shows  the  results  of  the  viscosity  calculation  for  both  rubidium 
and  cesium.  In  each  case  the  values  of  v  obtained  from  Eq.  (13),  as  well  as 
those  obtoined  from  Fontana's  calculations,  were  used  as  is  shown  on  the  graph. 
T  is  converted  to  degrees  Kelvin  as  explained  above.  Although  no  final  con¬ 
clusions  can  be  drawn,  the  evidence  seems  to  favor  the  larger  values  of  r. 

The  following  procedure  is  suggested  for  obUining  and  Cjj  when  C.j  of 
Eq.  (8)  is  known,  as  well  as  and  Calculate  €jj  from  Eq.  (11),  calculate 

C..  from  the  second  order  quantum  mechanical  perturbation  theory,  and  then 
JJ 

solve  for  (t..  . 

JJ 

It  may  appear  somewhat  strange  that  c  is  larger  for  rubidium  than  for 
cesium.  However,  there  is  apparently  no  theoretical  reason  why  this  should 
not  be  so.  These  resvilts  suggest  that  if  Rb2  and  Cs^  are  formed  by  "Liennard- 
Jones  binding,  "  Rb2  has  a  stronger  binding  energy  than  Cs2- 
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